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Abstract 

Measurements of the Higgs-boson production cross section at the LHC are an important 
tool for studying electroweak symmetry breaking at the quantum level, since the main 
production mechanism gg h is loop-suppressed in the Standard Model (SM). Higgs 
production in extra-dimensional extensions of the SM is sensitive to the Kaluza-Klein 
(KK) excitations of the quarks, which can be exchanged as virtual particles in the loop. 
In the context of the minimal Randall-Sundrum (RS) model with bulk fields and a 
brane-localized Higgs sector, we derive closed analytical expressions for the gluon-gluon 
fusion process, finding that the effect of the infinite tower of virtual KK states can be 
described in terms of a simple function of the fundamental (5D) Yukawa matrices. Given 
a specific RS model, this will allow one to easily constrain the parameter space, once 
a Higgs signal has been established. We explain that discrepancies between existing 
calculations of Higgs production in RS models are related to the non-commutativity of 
two limits: taking the number of KK states to infinity and removing the regulator on 
the Higgs-boson profile, which is required in an intermediate step to make the relevant 
overlap integrals well defined. Even though the one-loop gg ^ h amplitude is finite in RS 
scenarios with a brane-localized Higgs sector, it is important to introduce a consistent 
ultraviolet regulator in order to obtain the correct result. 



1 Introduction 



In the past decades there has been an enormous effort, both theoretically and experimentally, 
in trying to understand the origin of electroweak symmetry breaking. Within the Standard 
Model (SM), the Higgs mechanism provides a solution that demands the existence of a funda- 
mental self-interacting scalar field, the Higgs boson. After the analysis of significant amounts 
of data, the LHC and Tevatron experiments seem to observe, for the first time, an excess of 
events that could be associated with a Higgs-boson signal. If this excess is confirmed, a new- 
set of questions will need to be addressed. In particular, it will be crucial to understand how 
the electroweak scale M„cak is related, if at all, to other scales in nature such as the Planck 
scale Mpi. Searching for a framework that provides a natural connection between these scales 
has been one of the main reasons to propose theories beyond the SM. Some of these theories, 
such as supersymmetry, promote a perturbative extension of the SM up to scales close to Mpi. 
Others feature new strong interactions and provide a solution valid only up to energies a few 
orders of magnitude above the weak scale, without necessarily specifying the ultraviolet (UV) 
completion of the theory. 

In all cases, new particles are expected to appear at or slightly above the TeV scale to solve 
the naturalness problem of the Higgs mass parameter. These new particles are expected to 
have non-negligible effects on the production and decay rates of the Higgs boson. In the event 
of a Higgs discovery, any departure from SM expectations will then pave the way to probing 
models of new physics. From this perspective, precision Higgs-boson physics - much like rare 
weak decays and fiavor-changing neutral current (FCNC) processes - is an exquisite tool for 
probing the structure of electroweak interactions at the quantum level. The couplings of the 
Higgs boson to photons and gluons vanish at tree level in the SM, but they are non-zero at 
one-loop order and beyond. These couplings are therefore particularly sensitive to new heavy 
particles, which can propagate in the loops. 

In this article, we concentrate on analyzing Higgs-boson properties in models with a warped 
extra dimension, so-called Randall-Sundrum (RS) models [1], which represent a very appealing 
alternative to more traditional extensions of the SM, such as supersymmetry. These models 
provide not only a natural solution to the hierarchy problem, but also a very compelling 
theory of flavor. They feature a compact extra dimension with a non-factorizable anti-de 
Sitter (AdS) metric and two four-dimensional (4D) branes as the boundaries of the warped 
extra dimension. The AdS background generates an exponential hierarchy of energy scales, 
so that the fundamental mass scale near the "UV brane" is of order Mpi, while that near the 
"infrared (IR) brane" is suppressed by an exponential warp factor e~^'^'", with k being the AdS 
curvature and r the radius of the extra dimension. While all fundamental parameters are of 
Planck size, k ~ ~ Mpi, due to the curvature of the extra dimension the fundamental 
scale of the theory near the IR brane can be in the range of a few TeV. Hence, as long as 
the Higgs sector is localized near the IR brane, quantum corrections to the Higgs potential 
are cut off at the few- TeV scale and the hierarchy problem is solved. In addition, the large 
hierarchies observed in the spectrum of fermion masses and in the quark mixing matrix arise 
in a natural way from the localization of bulk fermions in the extra-dimensional space [2, 3, 4]. 
An important side effect of this mechanism is that flavor-changing couphngs between SM 
fermions, which arise due to the tree-level exchange of Kaluza-Klein (KK) gauge bosons, are 
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suppressed by the same small overlap integrals that generate the fermion masses [3] , a scheme 
referred to as the RS-GIM mechanism [5, 6]. 

In order to fulfill the constraints coming from electroweak precision data and FCNC tran- 
sitions, the low-lying KK excitations of the SM particles must have masses in the range of 
a few TeV, and hence may be out of the reach for direct production at the LHC. In recent 
years, many models have been proposed to alleviate these constraints and address the little 
hierarchy problem of the basic RS model with bulk fermions. The model building becomes 
rather involved and ranges from considering the effect of brane-localized kinetic terms for the 
gauge fields [7, 8] to enlarging the bulk gauge symmetries [9, 10, 11], thereby adding a whole 
new layer of additional matter fields. Another class of models takes the Higgs sector off the 
IR brane. This is particularly well motivated if the Higgs boson is identified with the 5**^ 
component of a gauge field ( "gauge-Higgs unification") [12, 13]. 

In this paper, we will concentrate on studying the effects of KK excitations on the pro- 
duction of a SM-like Higgs boson in the context of the simplest RS model containing bulk 
fermions and a minimal scalar sector localized on the IR brane. ^ In particular, we will scru- 
tinize the gluon-gluon fusion process, for which the presence of KK modes in loop corrections 
can have important effects. The type of considerations discussed in this work will also be rele- 
vant for the Higgs-boson coupling to two photons, which is the other loop-induced coupling of 
the Higgs field to gauge bosons. At present, the hterature contains two sets of results for the 
computation of Higgs production and decays in RS-type models that are at variance with each 
other. While in [14] the authors find a significant suppression of the gluon-gluon fusion cross 
section relative to the SM, the authors of [15] find an enhancement. Indeed, when applied to 
the same model, the two groups obtain effects on the hgg amplitude of approximately equal 
strength but opposite sign. In this article we revisit in detail both calculations and present 
alternative ways to derive them. For the first time, we obtain closed analytical expressions 
for both results in terms of the fundamental parameters of the model, valid to all orders in 
the ratio of the Higgs vacuum expectation value v and the typical mass scale Mkk of KK 
modes. We find that the calculations of [14, 15] are both correct from a mathematical point 
of view. The origin of their difference is related to a subtlety connected to the fact that, in 
order to compute its overlap with the fermion wave functions, it is necessary to regularize 
the Higgs profile in an intermediate step by taking it (slightly) off the IR boundary. The 
difference between the results obtained by the two groups arises from the different orders in 
which the Higgs regulator is removed and the sum over the KK tower is performed. We then 
show that the enhancement in the gluon-gluon fusion cross section obtained in [15] arises from 
very heavy KK modes. When the model is defined with a (warped) Planck-scale cutoff so as 
to solve the hierarchy problem, the cross section is suppressed, in accordance with the result 
obtained in [14]. 

The paper is organized as follows. In Section 2 we define our notation and present the 
setup of the problem. In Section 3 we present the low-energy Lagrangian for the effective hgg 
couplings induced by the virtual effects of KK fermions, which is valid for energies below the 
KK mass scale. We also describe the setup of the two calculations, one summing first over the 

^Extensions to models with a custodial symmetry in the bulk, or generalized warped extra-dimensional 
models wih be considered elsewhere. 
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infinite KK tower and then taking the Higgs regulator to zero, and the other considering the 
hmits in the reversed order. In Section 4 we perform the first calculation and reproduce the 
results of [15] in a different way by using the five-dimensional (5D) fermion propagator of the 
theory. In this way, we generalize the result to all orders in v/Mkk- In Section 5 we perform 
the second calculation, considering first a simplified case, in which it is straightforward to 
obtain a closed answer for the sum over KK modes. Based on symmetry considerations, we 
then make a conjecture for how to extend this result to the general case, which reproduces 
the numerical results obtained in [14]. In Section 6 we explore in more detail the origin of the 
difference between the two approaches, which in RS models with an IR-localized Higgs sector 
arises from unphysical contributions of KK modes with "trans-Planckian" masses. We then 
argue in Section 7 that the use of a consistent UV regularization scheme, which for instance 
is provided by the physical UV cutoff inherent in RS models, eliminates these contributions, 
thus favoring the result corresponding to a suppression of the gluon-gluon fusion cross section. 
Phenomenological imphcations of our findings for Higgs-boson searches at the LHC are studied 
in Section 8, while Section 9 contains our conclusions. 

2 Preliminaries 

We work with the non-factorizable RS geometry [1] 

ds^ = e-'^^"^) r^^, dx^dx" - ; a{(j)) = A;r|0| , (1) 

where denote the coordinates on the 4D hyper-surfaces of constant with metric rj^iy = 
diag(l, —1, —1, —1). The 5*^ dimension is an /Z2 orbifold of size r, labeled by G [— tt, tt]. 
Two 3-branes are located at the orbifold fixed points at = (UV brane) and ^ = tt (IR 
brane). The curvature parameter k and the radius r of the extra dimension are assumed to 
be of Planck size. The warped Planck scale Mpi e"""^*^-* sets the effective, position-dependent 
fundamental scale at a given point along the extra dimension. It serves as a natural UV 
cutoff, since quantum gravity would become relevant above this scale. In order to solve the 
(big) hierarchy problem, one assumes that the effective cutoff scale on the IR brane, Atbv = 
Mpi e~^*^'^^ = Mpi e, is in the range of 20 TeV or so. This requires that L = km = — \\ie 34. 
The "little hierarchy problem" , the fact that the Higgs-boson mass is two orders of magnitude 
smaller than the cutoff scale, is not addressed in the minimal RS framework. The warped 
curvature scale, Mkk = ke~"'^^^ — ke, is assumed to he somewhat lower, in the range of a few 
TeV. It sets the mass scale for the low-lying KK excitations of the SM fields. For instance, 
the masses of the first KK photon and gluon states are approximately 2.45 Mkk [16]. It will 
be convenient to introduce a coordinate t = ee'^^'^\ which equals e ~ 10~^^ on the UV brane 
and 1 on the IR brane [2]. 

We consider the minimal RS model in which all SM matter and gauge fields propagate 
in the bulk, while the Higgs boson hves on the IR brane. The 5D theory contains three 
generations of massive fermions Q = {U,D)'^ and u, d, which transform as doublets and 
singlets under the SM gauge group SU{2)i, respectively. The elcctrowcak gauge symmetry 
SU {2)lX U (l)y is broken by the Higgs sector on the IR brane. The 5D fields are functions of 



3 



X and 0. The doublet fields Q have left-handed zero modes, while the singlet fields ti, d have 
right-handed ones. These correspond to the chiral fermions of the SM. 

In this paper we are particularly interested in the couplings of the SM Higgs boson to 
fermions. In unitary gauge, the Higgs-boson couplings to SM quarks and their KK excitations 
read 

A,,--EE^^n^?-f^??^ + h.c., (2) 

q=u,d m,n 

where the Yukawa couplings are given in terms of the overlap integrals [14] 

(3) 

and likewise in the down- type quark sector. Here n labels the different mass eigenstates q^^\x) 
in the 4D effective theory, such that n = 1,2,3 refer to the SM up-type quarks u,c,t, while 
n — 4, . . . , 9 label the six fermion modes of the first KK level, and so on. The quantities 
Ci^\t) and Si'^\t) with A = Q,u,d are diagonal 3x3 matrices in fiavor space, which 
contain the Z2-evcn and odd fermion profiles along the extra dimension, respectively. These 
can be expressed in terms of combinations of Bessel functions, whose rank depends on the 
bulk mass parameters Cq = Mg/k and Cu^a = —Mu^^/k of the 5D fermion fields [2, 3]. The 
SU (2) L gauge symmetry in the bulk implies that SU (2)-doublet quark fields have common 
Cq parameters. The 3-component vectors on the other hand, describe the fiavor mixings 
of the 5D interaction eigenstates into the 4D mass eigenstates, which are generated by the 
Yukawa interactions on the IR brane [14]. Because of electroweak symmetry breaking, these 
vectors are different for A = U,D^ u, d. For simplicity, from now on we use the generic notation 
Q for U, D and q for u, d. The 3x3 matrices Yq contain the dimensionless Yukawa matrices 
of the 5D theory, which contrary to the SM are assumed to have an anarchical structure, i.e., 
they are non-hierarchical matrices with 0{1) complex elements. The hierarchies of the Yukawa 
matrices of the SM quarks in the effective 4D theory are explained in terms of a geometrical 
realization of the Froggatt-Nielsen mechanism in RS models [17, 18, 19, 20]. 

In the second line of (3) we use a compact notation, in which the profile functions of the 
left-handed (right-handed) interaction eigenstates that can mix into the left-handed (right- 
handed) components of the 4D mass eigenstates are collected in 6-component vectors (t) . 
These vectors obey the orthonormality conditions (with A — L,R) 

fdtQf^\t)Qf{t)^5mn. (4) 
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They are defined in terms of tlie KK decompositions 



/-\(n)/-,\ (n / \ 



(5) 



where q^"]i(x) = | (1 ^ 75) q^"'\x) denote the left- and right-handed components of the n*^ 
Dirac fermion. 

Note that even the diagonal Yukawa couplings g"^^ derived from (3) are in general complex 
numbers, such that pseudo-scalar currents appear in the Lagrangian 

^tl = - E E [^<9L) h + Im(^L) h i75 g^"^] , (6) 

q=u,d n 

which in principle could have interesting implications for phenomenology. Unfortunately, 
however, for the SM fermions the imaginary parts of the Yukawa couplings are very strongly 
suppressed [14]. 

It has been emphasized in [21] that in order to properly evaluate the Yukawa couplings it 
is necessary to regularize the Higgs-boson profile. The reason is that with a brane-localized 
Higgs sector the odd fermion profiles Sii^\t) are discontinuous at t = 1, and hence the overlap 
integral of a product of two such functions with the naive Higgs profile 6{t — 1) is ill defined. 
In (3), wc have therefore replaced this profile with a regularized delta function 6'^{t — 1). The 
precise shape of this function will be irrelevant; however, it is important that the Higgs profile 
is non-zero only within a small interval of width 77 <^ 1 next to the IR brane, 6'^{t — 1) 7^ 
only a 1 —T) < t < 1, and that it has unit area. The r^-dependence of the Yukawa couplings is 
implicit in our notation, but it will play an important role in our analysis. 

In the presence of the regularized Higgs profile, the bulk equations of motion (EOMs) for 
the profile functions read [14] 

I (7) 
-| Qt\t) = -xn Qf{t) + M,{t) Qt\t) , 
where Xn = mg^/MKK are the mass eigenvalues, and 

1 /cq \ V ^, / 



is the generalized mass matrix. Without loss of generality, the hermitian bulk mass matrices 
Ca of the 5D theory can be taken to be diagonal. The boundary conditions are such that the 
odd profiles vanish on the two branes, which implies 

(0 l)Qj^)(t,) = 0, (1 O)qP{U) = 0; for t, = {e,l}. (9) 
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Note that these simple Dirichlet boundary conditions only hold because the Higgs profile has 
been regularized and is no longer singular on the IR brane. When t is lowered away from 1, the 
profile functions change rapidly over the small interval in which the Higgs profile is non-zero. 
A careful analysis of the solutions in the vicinity of the IR brane shows that at i = 1 — 77 the 
solutions obey the equation (assuming 77 <S 1) [14] 

(vsk^' i)eS;"(i-.)^o. (1 -^n)eS'(i-.)-o. (10) 



where 



is a positive definite, hermitian matrix given in terms of the 5D anarchic Yukawa matrices of 
the model. Note that the mixed boundary conditions (10) hold irrespective of the shape of 
the regularized Higgs profile, as long as it is confined to a box of width rj. It is thus possible 
to take the limit ^ corresponding to a brane-localizcd Higgs boson, drop the Yukawa 
couphngs in the bulk EOMs (7) and (8), and impose the mixed boundary conditions on the 
IR brane, i.e., aX t — 1~, where 1~ denotes a point infinitesimally to the left of the IR brane. ^ 
Note, however, that to derive the result (10) it was important to regularize the Higgs profile 
in an intermediate step [21]. 

The mixed boundary conditions (10) evaluated at t = 1" are sufficient to determine the 
masses and profiles of the fermionic KK modes required, e.g., for the analysis of many tree- 
level flavor-changing processes [19]. However, when computing overlap integrals of fermion 
profiles with the wave function of the Higgs boson the correct procedure is to start from a 
regularized Higgs profile, compute the relevant overlap integrals, and then take 7] to zero [21]. 
Using the explicit expressions for the profile functions near the IR brane derived in [14] , which 
were obtained from the solution of the system of equations (7) to (9) assuming 77 <^ 1, it is 
straightforward to derive from (3) that 



9ln - ^ [dt 5\t - 1) aL^)t C(?)(l - 77) ^ 

cosh^(^''(i - 1) Xq) - sinh^(^^(i - 1) Xq) 



(12) 



where the cosh and sinh functions in brackets arise from the Z2-even and odd fermion profiles 
in (3), respectively. We have introduced the notation 

^'?(t-l) = [ (is5'?(s-l) (13) 



for the integral of the regularized Higgs profile, such that O'^iQi) = and 9^{t — 1) = 1 for 
t < 1 — rj. We observe that the products of the i-dependent fermion profiles add up to a 



^For rj ^ the profile functions are discontinuous on the IR brane, so it is necessary to stay an infinitesimal 
amount away from the boundary. 
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Figure 1: Effective hgg couplings induced by the exchange of virtual KK quarks. 



constant, so we are left with the normalization integral over the Higgs-boson wave function, 
which equals 1. We thus obtain 



cosh X, 



2^ Y.Ci^Hl )a 



(14) 



Interestingly (and somewhat unexpectedly), the same result can be obtained in a more 
naive way. By evaluating relation (3) with a brane-localized Higgs boson, we obtain 



lim gl 



TT 



Le 
V27r 



,(Q)t 



iQ) 



,(Q)t 



(15) 



where in the second step we have used the mixed boundary conditions (10) at t = 1^. This 
result agrees with (14) by virtue of the definition of Yg in (11). 

The Yukawa couplings derived in (14) will be used in the analysis in Section 5. For the 
discussion in other parts of our paper, it will be necessary to keep the regulator rj non-zero 
until the sum over the tower of KK modes has been performed. 



3 Effective low-energy theory for hgg couplings 

We are now ready to derive the effective low-energy Lagrangian for the Higgs-boson couplings 
to a pair of gluons, which are induced by the exchange of KK quarks. This Lagrangian is valid 
at energies below the scale Mkk, at which these states can be integrated out. The relevant 
Feynman diagram arising at one- loop order is shown in Figure 1. Since the gluon couplings to 
fermions are diagonal in the mass basis, a single quark state g'-"-' runs in the loop. Summing 
over the KK tower, we obtain 

^hgg = Ci (/x) -j^ - G^, G^' -C, ifi) - G^, G^ , (16) 

where ^ < Mkk is the scale at which the effective operators are renormalized, and C^^'" = 
_l^^^lya|3 Qo^^ (with e°^^^ = —1) is the dual field-strength tensor. Throughout this paper, v 
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denotes the value of the Higgs vacuum expectation value in the RS model, which differs from 
the SM value vqm ~ 246 GeV by a small amount [22]. To first order in v'^/M^^, we obtain 
from the shift of the 1^-boson mass predicted by the RS model [19] 

"'^ ^L-l + ^), (17) 



where L f» 34. At the matching scale fj, — Mkk, we find for the Wilson coefficients entering 
(16) the one-loop results 



q=u,d n ^" 

cr(MKK)=EE'^^5(-..). 



q=uA n ^" 

(18) 



q=u,d n ^" 



where Tq = 4to^ /ml — is. The sums extend over all KK fermion states but exclude the SM 
fermions {n ^ 1,2,3). This is indicated by the prime on the sum symbol. Using the results 
from [23, 24], we arrive at 

Air) = — fl + (1 - r) arctan^ , ^ 1 , Bir) = r arctan^ , ^ . (19) 

Since the KK quarks are much heavier than the Higgs boson, it is sufficient to use the asymp- 
totic values A{oo) — B{oo) — 1 of these functions. It then follows that 



Cf ^(Mkk) + i C5^^(Mkk) = E E — ■ (20) 

q=u,d n ^" 
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The real part of the sum determines Cf^, while the imaginary part gives C^^. 

We recall at this stage that each term in the sum (both the Yukawa couplings and the 
mass eigenvalues) depends on the regulator rj used to smear out the Higgs profile. In addition, 
one needs to worry about the convergence of the infinite sum over KK modes. Indeed, since 
the Yukawa couplings g^^ are of 0{1) and the masses of the KK modes are approximately 
evenly spaced multiples of the KK scale, naive dimensional analysis would suggest that the 
sum diverges logarithmically. In order to define it properly, one should therefore regularize 
the sum, for instance by introducing a cutoff on the highest KK level that is included. We 
thus define 

3+eN q 

^«(^'^)= E (21) 

where the sum now includes the SM quarks (n = 1, 2, 3) plus the first levels of KK modes. 
Each KK level contains six modes with 0{v) mass splittings, while the different levels are spht 
by an amount of O(Mkk)- It will turn out that the sum over modes is finite despite naive 
expectation, because non-trivial cancellations happen among the six modes contained in each 
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KK level. Nevertheless, the notion of a cutoff on the KK level will play an important role in 
our analysis. 

The Wilson coefficients in (20) are now obtained as 



q=u,d L n=l,2,3 * 



(22) 



where the last term subtracts the contributions from the SM quarks, each of which equals 1 up 
to higher-order corrections in f/MxK- These contributions can be derived using the explicit 
expressions for the profile functions of SM quark fields near the IR brane given in [14] as well 
as the relations for the aif -' vectors presented in [19]. We obtain 



E 



V9nn 



Tr 



n=l,2,3 



2X. 



sinh 2X. 



— e 



1 ' 



(23) 



where the quantity 



= Tr {bQ + ^,) + O 



33 



(24) 



contains some small corrections to the leading term. It receives its main contributions from 
the third-generation terms, because the quantities (5(5,g)„„ oc defined in [19] are 

chirally suppressed for all light quarks. To a very good approximation 



33 



1 A 1 



3i\ 



33 



2a 



U3 



3 + 2c„. 



(25) 



and (5ji)33 is given by the same formula with c^^ — )■ CQ^ and (1^)3^ — >■ {Yv)i3- Analogous 
expressions hold in the down-type quark sector. The function 



F\c) 



1 + 2c 
1 - ei+2c 



(26) 



determines the values of the fermion profiles on the IR brane [2, 3]. Note that in the approx- 
imation used here both Eq and the zero-mode sum in (23) are real and do not contribute to 
the Wilson coefficient Cf^ in (22). 

The cross section for Higgs-boson production in gluon-gluon fusion can now be written as 



(^{gg ^)rs 



I I _L I I 



(27) 



where Ky has been defined in (17), while 
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K>5g 



(28) 



Here K,t = vg^-^/rrit and = vg^^/rrib encode the modifications of the Higgs-boson couphngs 
to top and bottom quarks with respect to their SM values. The former quantity is given to 
high accuracy by [25] 

7)2 (YY^'Y) 

^^^^-^ ]y\ "^'^ - ('^^)33 - iSu)s, , (29) 

and an analogous expression with u ^ d and U ^ D holds for k^. 

The Wilson coefficients must be evolved from the high KK scale to the Higgs-boson 
mass scale. Their scale dependence is governed by evolution equations. If the KK contributions 
are dominated by the corrections due to the lowest-lying modes, it is practical to integrate out 
all KK modes at a common scale of order Mkk, as already done in (18) above. The evolution 
from this scale down to the weak scale is described by the rcnormalization-group equations of 
the SM. For the QCD evolution of the Wilson coefficient Cf ^(/x), one obtains the exact result 
[26, 27] 

gf^(/^) ^ I3{as{^))/al{^) ^13 a,(/.) - ^.(Mkk) . . 

CfHM^^) /3{asiMKK))/c^{MKK) 14 tt ■■■■ ^ ^ 

In practice, the evolution from the KK scale of several TeV down to a; ruh has only a small 
effect of about 1% on the value of Cf^, since no leading logarithms appear in this result. 
Since the operator multiplying the coefficient Cf^ in (16) is connected by the Adler-Bell- 
Jackiw anomaly to current operators with vanishing QCD anomalous dimension, it follows 
that this coefficient is scale independent in QCD, i.e., C|^^(/i) = C|^^(Mkk)- 

Much of the analysis in the present paper is concerned with the fact that the order in 
which the limits rj ^ and A^" — >■ oo in (22) are taken is not irrelevant, because these two 
limits do not commute. In our previous work [14], we have evaluated the Yukawa couplings 
in (3) mode by mode, taking the regulator t] to zero after computing the relevant overlap 
integrals. We have then numerically evaluated the contributions of the first few KK levels 
to the Higgs-boson production cross section, observing that the sum over modes converges 
and approaches a limiting value after the summation over several KK levels. In essence, this 
approach corresponds to taking the limit 77 — > first, and hence 



3+6Af 

WCGHNP)^ hm [limE,(7V,r^)l = lim V ''^^ 



n=l 1^ 



(31) 



Yet we did not attempt to derive an analytical expression for the infinite sum. In the approach 
taken in [15], on the other hand, one considers the infinite sum over modes from the very 
beginning. This is accomplished by means of completeness relations for the fermion profiles. 
The regularization of the Higgs profile is taken to zero at the end of the calculation. It follows 
that in this approach 

00 „ 



.N^oo 



Ei^^^) = lim lim EJN,ri) = lim V ^ . (32) 



^ mq„ 



In the following two sections, we will derive closed analytical expressions for the limiting values 
in (31) and (32). In Section 6, we will then discuss the physical interpretation of these results. 
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4 Summing first over the infinite KK tower 



We begin with the study of the sum in (32), in which the regulator on the Higgs profile 
is removed after the infinite sum over KK modes has been performed. The authors of [15] 
succeeded to evaluate this sum using the completeness relations (^4 = L, R) 

J2QA\t)Q^f{t')^5{t-t') (33) 



for the fermion profiles. The analysis relied on a perturbative treatment of the Yukawa cou- 
plings, which yields the first non-trivial contribution in the expansion in powers of v/Mkk- 
A subtlety in this calculation, worth recalling, is that in the perturbative approach the odd 
fermion profiles S^\t) vanish on the IR brane. As a result, the contribution of each individ- 
ual mode, and hence of any truncated KK sum, vanishes in the limit where the width of the 
Higgs-boson profile is taken to zero. A non-zero result for the (9(f^/M^j^) correction to (32) 
is obtained only after summing over the infinite tower of KK states. In our notation, the main 
result of [15] reads 

ElATZ)^Xr(^l + 5 + ...V (34) 



An elegant way to derive a closed expression for the sum (32), valid to all orders in v/M^k, 
is to relate it to the 5D fermion propagator of the theory. We study propagator functions built 
using the 6-component spinors 



(35) 



where the profile functions Q^^\t) have been defined in (5). The grand, 6x6 propagator in 
the mixed momentum/position representation [28, 29] is 



{t,t';p)^ J d'xe'P'' {0\T{QLit,x) + QR{t,x)){QL{t',0) + QR{t',0)) |0) 



E 



m. 



(36) 



(37) 



where T denotes time ordering. In our notation, the Dirac operator takes the form 

d 

P = # - Mkk 75 ^ - Mkk Mq{t) , 

with the generalized mass matrix defined in (8). Using the bulk EOMs (7) and the complete- 
ness relations (33), it is straightforward to show that 

VS'>{t,t']p) = S{t-t'). (38) 

Since with our exact treatment of Yukawa interactions there are no massless zero modes of 
the fermion fields, the propagator does not exhibit a singularity at = 0. We can therefore 
study the special limit p^ ^ without complications. We obtain 



S'^{t,t';0)^- 



(39) 
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where 



(40) 



It now follows from the second relation in (3) that the infinite sum in (32), for fixed non-zero 
77, can be expressed as [30] 

"1 



where 



lim E,(7V,77) 

AT— >-oo 



^lL(^,0 = ^Tr 



dt5\t-l)Tl^{t,t), 









(41) 



(42) 



Notice that only the off-diagonal blocks of the propagator enter in this result. 

The Dirac equation (38) for the 5D propagator implies the differential equation 



d_ 
di 



M,{t) 



1 



Ml 



KK 



5{t - 1') . 



(43) 



We will construct solutions to this equation assuming first that t distinguishing the cases 
where t > t' and t < t'. Later, these solutions are patched together by means of the jump 
condition 

'AUt + d,t)-A%^{t-d,t)\ ^ 



lim 



Ml 



KK 



(44) 



Since we have regularized the Higgs profile, the boundary conditions on the UV and IR branes 
take the simple form 

(1 0)AUe,O = (l 0)AU1,0 = (0 0), (45) 

which follows since the Z2-odd profile functions S^\t), which sit in the upper components of 
the right-handed spinors in (5), vanish on the branes. The three equations above determine 
the matrix function A|j^(t, t') completely. 

For our purposes, it suffices to construct the solution in the region close to the IR brane, 
where I — rj < t,t' < 1. For t > t', equation (43) is solved by the t-ordered exponential 

ft 



\t>t' 



Texp 



^ dsMq{s) 



(46) 



where (45) implies that the boundary function A|j^(l,i') has only lower components. The 
ordering symbol is required since the matrices 



dsMq{s) = Int 



Cq 
-c„ 



V2M] 



^"(t-l) 



KK 



Yq 





(47) 



with 9'^{t — 1) as defined in (13), do not commute at different t values. For coordinates t near 
the IR brane, however, the first term in (47) can be neglected, and we can approximate 



y2M] 



9'\t 



KK 



Yq 

i;t 



+ 0{ri) 



(48) 
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where the harmless 0{ri) terms can be dropped. The resulting matrix exponential can then 
be evaluated in terms of hyperbolic trigonometric functions. We then lower the value of t, 
implement the jump condition (44) when t crosses t', and finally derive an expression for 
A|j^(l — r),t'), which still contains two unknown coefficients that are matrices in generation 
space and functions of t'. They can be determined using the boundary condition on the 
UV brane. For t < 1 — rj the generalized mass matrix Aiq{t) in (8) is diagonal, since the 
contribution from the Higgs profile vanishes. The first boundary condition in (45) therefore 
imphes that (1 0) A|j^(l — r/, t') = (0 0) as long as t' > 1 — rj, and this equation yields the 
desired constraints on the propagator function. 

In (42) we only need the off-diagonal blocks of the function A^^(i, t') for 1 — 77 < f < 1. 
We find 



rlJi,0 = Tr 



sinh X„ 



cosh ( Xn 1 



ds S^is-l) 



(49) 



where i< = min(t,t') and t> = max(t, t'). Notice that, for the special case where t — t' , the 
result in expression (49) is independent of both t and 77. This means that in (41) only the 
normalization of the Higgs profile enters, which is independent of its shape and equal to 1. 
We conclude that the sum in (41) is in fact independent of 77, and hence 



y;(ATZ) 



Tr(XgCothXg) = Tr 1 + 



45 



(50) 



The first non-trivial term in the Taylor expansion agrees with (34). Note that this result is 
real, and hence it only contributes to the Wilson coefficient Cf^. Remarkably, the infinite 

sum of KK states converges despite of the fact that it is superficially divergent. This hints to 
the existence of intricate cancellations between different contributions to the sum. From (22) 
and (23), we now obtain 



Cf ^(Mkk) = ^ [Tt {X, tanh X,) + Re £,] 



q=u,d 

E 

q=u,d 



X^ 

Tr(X2-^ + ... ) + Reeg 



(51) 



and 



C5™(Mkk) = Im(£„ + £d) 



(52) 

Recall that the correction terms Eq are very small, and to an excellent approximation they are 
real, see (24). The Wilson coefficient Cf^ is therefore predicted to be tiny, of order v'^/M^j^ at 
most. The contribution to the coefficient is positive and hence yields to an enhancement 
of the gg ^ h amplitude compared with the SM. 



5 Removing first the regulator on the Higgs profile 

We now consider the sum in (31), for which compared to (32) the limits N ^ 00 and ?] — )■ 
are taken in the opposite order. In this case one first considers a finite sum over KK modes 
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and derives the relevant Yukawa coupling for each mode by regularizing the Higgs profile, 
computing the overlap integral in (3), and then taking the limit — )■ 0. The relevant Yukawa 
couplings in that limit have been given in (14). 

Because the sum in (31) extends over a finite number of KK levels, it is not possible to use 
the elegant method of 5D propagators described in the previous section. In order to obtain 
a closed expression nevertheless, we adopt the strategy of first finding a solution in a special 
case, where the bulk EOMs and the eigenvalue equation determining the masses of the KK 
modes can be solved analytically. We will then argue that, as in the case (50), the solution is 
independent of the bulk mass parameters. 

The special case we consider is that of one generation of fcrmion fields, whose bulk mass 
parameters vanish: cg = c„ = q = 0. The EOMs for the profile functions are then solved 
in terms of simple trigonometric functions, which can be evaluated in the limit e ^ (recall 
that e 10"^^ is tiny in the RS model). We find that the KK masses in units of Mkk are 
given by the solutions to the eigenvalue equation 



tan^Xn = tanh X, 



1 1 



V2M] 



(53) 



KK 



where without loss of generality we assume that Yq is real and positive. This equation can be 
solved to give 



n — 1 



TT + xi ; n = 1, 3, 5, 



n 



(54) 



TT — xi ; n = 2, 4, 6, 



where xi — arctan(tanh Xq) denotes the mass of the zero mode (the "SM quark" ) in units of 
the KK scale. The corresponding even profile functions are 




CP(i)aW) = cos(x„t), 




C(^)(t)ai^) = ±cos(x„i), 



(55) 



where the upper (lower) signs hold for odd (even) values of n. Inserting these results into (14) 
and using (54), it follows that 



i>nn 



±1 Y„ 



±1 



cos Xr. 



(56) 



^/2 cosli^Xq " ^/2 cosh 2Xg ' 

We can now readily compute the sum over KK modes required in (31), with the result that 

i+2Ar 



\imEq{N,ri) 

rj—^O 



1 gcncr. 
CA=0 



= E 

n=l 



Xn 



cosh 2X„ 



1 ^ 



1 



kn + xi kn — xi 



(57) 



In the one-generation case, n = 1 refers to the zero mode, while n >2 labels the KK modes. 
The first term in the bracket of the final expression arises from the "SM quark" , while the sum 
is over pairs of KK modes belonging to the A;*'^ KK level. Note that for large k the individual 
terms in the sum fall off only like 1/k, but each pair combines to a contribution decreasing 
hke l/k"^. Hence, as long as we sum over complete levels of KK states first, the sum over k is 



14 



convergent, and it is possible to take the limit N — )■ oo. We thus obtain for (31) in the case 
of a single fermion generation with vanishing bulk masses 



j^(CGHNP) 



1 gener. X coth X 



^ CA=o cosh 2Xq 

which differs from the finite sum in (57) by terms of 0{1/N). 

For one generation, it is not difficult to solve the EOMs numerically also for general ca 7^ 
0. The profile functions and eigenvalue equation in this case involve Bessel functions. We 
have confirmed in this way that the simple answer (58) still holds in this more general case. 
Given the similarity of the above result with equation (50) derived in the previous section, we 
conjecture that in the case of three generations 

,,ec™.,^,(W).,(,_^,H^,...). 

Once again the answer is real, such that only Cf^ receives a contribution. Combining (22), 
(23), and (59), we then find 



q=u,d 

= E 

q=u,d 



^ , XotanhXo , ^ 
-TV( ^ ^ 1+Re£, 
cosh2A„ 



^1 



Tr ( + ^ + ) + Reea 



(60) 



while is still given by (52). This time the contribution to is negative and thus yields 
a suppression of the gg ^ h amphtude compared with the SM. We have checked that this 
formula indeed reproduces our numerical results obtained in [14]. Incidentally, the first term 
in the expansion in powers of has the opposite sign from the result (51), which was first 
obtained in [15]. 

The conjecture that equation (58) can be generalized to the case of three generations will 
be proved elsewhere. We note, however, that symmetry arguments help us to constrain the 
form of the answer. It is plausible to assume that, after the limit N ^ 00 has been taken, the 
result for the sum (31) should depend only on the fundamental parameters of the underlying 
5D theory, i.e., the bulk mass parameters and Yukawa couplings. All reference to the properties 
of individual states, such as their masses or profile functions, should disappear. We express 
this assertion by writing 5]^^*^^'^^^ = 'E(Yq,CQ,Cq). The dependence of this quantity on the 
fundamental parameters can be constrained using symmetry arguments. The EOMs (7) and 
boundary conditions (9) are valid in an arbitrary basis, in which the bulk mass matrices Ca are 
not necessarily diagonal. These relations are invariant under a set of three global symmetries. 
The first one is a symmetry under the exchange of SU{2) doublets and singlets along with 
Yq •<-> . In addition, there are two symmetries related to unitary transformations of the 
Yukawa and bulk mass matrices. When combined with the fact that in the one-generation 
case the result is found to be independent of the bulk mass parameters, these symmetries 
imply that the quadratic term in Xq must indeed be of the simple form shown in (59). Using 
symmetry considerations alone, we can however not exclude a dependence on the parameters 
Ca starting at C(X^), provided that it cancels in the case of one generation. 
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6 Reconciling the results 



The discussion of the previous two sections shows that the calculations presented in [14] 
and [15] are free of mistakes. The results for the effective hgg couplings derived in these 
two papers differ, in magnitude and sign, because both groups considered different orders in 
which the limits ?7 ^ and ^ oo in (22) were taken, and the two limits do not commute. 
Hence, the question poses itself which of the orders of hmits is the more reasonable one from 
a physical point of view. We will address this question in Section 7. 

A second, equally puzzling question arises once we realize that the non-commutativity of 
the limits ?7 — ?■ and N ^ oo implies that the origin of the discrepancy between the results (50) 
and (59) must be due to contributions from KK modes with large n and hence very heavy 
masses, and related to the behavior of the profile functions for such modes very close to the 
IR brane. The fact that in both approaches one finds that the infinite sum over KK modes 
converges, signaling that heavy modes decouple, raises the question: how is it possible that 
modes with very heavy masses give an 0{1) contribution to the sum, which is large enough 
to change the sign of the answer? 

In order to understand the physics behind this effect, it is useful to study a toy model 
that can be solved exactly also for finite rj. To this end, we consider again the case of a 
single fermion generation and vanishing bulk mass parameters, which proved so useful for our 
analysis in Section 5. To obtain analytic solutions to the bulk EOMs for the fermion profiles 
we adopt a particularly simple form for the regularized Higgs profile, namely a box of width 
77: ^ 

{-; for l-?7<t<l, 
V (61) 
; for e < i < 1 - ?7 . 

Moreover, since for ca — the wave functions are non-singular near the UV brane, it is 
possible to set e = 0. In Appendix A we derive the explicit forms of the bulk profiles for 
this toy model and the eigenvalue equation that determines the masses of the KK modes. 
We then present a formula for the sum Eg(iV, ?]) entering (22). It turns out that the nature 
of the solution differs depending on whether a;„ < z 01 Xn > z, where we have introduced 
the abbreviation z — vYg/ {y/2MKKV) — ^g/v- The latter condition corresponds to masses 
niq^ > vYq/{y/2rj) — Mweak/^, whcrc Mweak = vYqj \f2 is of order the weak scale, since we 
as usual assume Yq = 0{1) in the anarchic RS model (recall that we can choose Yq real and 
positive). The appearance of the scale M^^ak/v^ which for very small rj lies far above the TeV 
scale, will be of crucial importance to the resolution of the puzzle. 

We find that the mass eigenvalues in our toy model are determined by the conditions 



tan [xn{l - rj)] ^ ^ 



±\l ^ ^ tanh [y/ — x\ ry] ; for x„ < 2; , 

Z 31 X^i 

(62) 



Xn i Z 



tan [y^x'l - z'^ 77] ; for x„ > 2; , 



where as before the upper (lower) signs hold for odd (even) values of n. These equations can 
easily be solved numerically. The sum over the KK contributions to the effective hgg couplings 
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is then obtained by evaluating the sum 

1+2N „ 1+2N r 



„_i '"'Qn „_i L \ -^n / 



n=l n=l 

= 1 + 



(63) 



±^sin [2x„(l -77)] -r]{zTXnCOS [2x„(l -ry)]) 



In the left plot in Figure 2, we show numerical results for the sum (63) as a function of the 
number of KK states included, for four different values of 77. For the purpose of illustration 
we take Xq — 0.5. We observe that for low values of N the sum quickly converges toward a 
value close to the result (31), which equals 0.701 in the present case. But then there is an 
intermediate region, roughly given by the range O.l/r] < N < 10/r], in which the value of 
the sum changes by an 0{1) amount. After this transition region, the sum converges to the 
value corresponding to the result (50), which equals 1.082 for our choice of Xg. In order to 
understand the origin of the three regions - in particular the appearance of the intermediate 
region, in which the sum grows by an 0(1) amount despite of the fact that the corresponding 
KK masses are extremely heavy - we show in the right plot the values of the Yukawa couplings 
5f^„ in our toy model (for the case where t] = 10"'^). For low values of n, the couplings for each 
pair of KK states in the same KK level have equal values and opposite signs. In this case, like 
in (57), the contribution from the k^^ KK level is 



Xg f 1 1 \ _ Xg 2xi 



cosh 2Xq ykn + xi kn — xi J cosh 2Xq k^n'^ — x 



(64) 



In this region the sum over k is convergent and dominated by the contributions from the 
first few KK levels. The situation changes drastically in the intermediate region, where the 
average coupling (y'avg(fc) in each KK level, shown by the dark points in the right plot, no longer 
vanishes. There is thus a range of logarithmic growth of the sum, over which its contribution 
can be estimated as 

y^^i^ln^, (65) 

where (gavg) denotes the mean value of gavg{k) in the interval k G [A^i,A^2], in which the 
average coupling in each KK level departs from zero. For even larger values, the average 
coupling decreases quickly, giving rise to very small contributions to the sum. 

Physically, the intermediate region arises because, for fixed rj, there exists a minimum 
KK mass rriq^ ~ M„cak/v beyond which the KK profile functions begin to penetrate the 
box modeling the Higgs-boson profile. When this happens, the cancellation of the Yukawa 
couplings of KK modes within one KK level is no longer operative. Only for yet much higher 
KK levels, when the profiles exhibit a large number of oscillations within the box, the couplings 
average out to zero and hence decrease with increasing k. While this behavior might seem 
strange at first sight, one should remember that by naive dimensional analysis the sum T,q{N, 77) 
is logarithmically divergent for large N. It is only due to subtle cancellations that it converges; 
however, when the profiles penetrate the box these cancellations no longer occur, and the 
generic logarithmic growth arises. 
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10 100 1000 10^ 

Figure 2: Left: Partial sums Sg(A^, rj) for different values of rimax = 1 + and 1] in the 
toy model with one generation and vanishing bulk mass parameters. The curves refer 
to ?7 = 10~^ (red), 10~^ (green), 10~^ (blue), and (black). Right: Yukawa couplings 
g'l^ of the KK fermion states in the toy model for t] = 10~^. The dark points show the 
Yukawa couplings averaged over the pair of modes in each KK level. 



7 Relevance of the UV regulator 

In order to decide which of the two calculations presented in [14] and [15] is physically more 
meaningful, we recall the importance of using a consistent UV regularization scheme when 
evaluating the gluon-gluon fusion amplitude. This is true in the SM, and even more so in its 5D 
extensions. Here we will study two different regularization schemes: dimensional regularization 
and the use of a hard momentum cutoff. While a dimensional regulator is particular convenient 
in that it preserves gauge and Lorentz invariance, the second option is also a natural choice 
in the present case. This is because the RS model must be considered as an effective theory 
below the Planck scale, which requires a UV completion incorporating the effects of quantum 
gravity. A peculiar feature of warped extra-dimension models is that the effective UV cutoff 
depends on where the theory is probed along the extra dimension [31] (see also [32, 33, 34, 35] 
for related works). The physical reason is that due to warping the fundamental length and 
energy scales change along the extra dimension. More specifically, the effective cutoff scale at 
the position t in the extra dimension is of the order of the warped Planck scale 



Auv(t) ~ Mpie 



A. 



TeV 



PI 



t 



t 



(66) 



The cutoff should be sufficiently high that at least a small number of KK modes have masses 
below Axev, and hence Atcv/Mkr = C(10) or so. Otherwise incalculable "threshold cor- 
rections" of order M^-^/A'^^y become important and take away the predictive power of the 
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model. We stress that imposing a UV cutoff is crucial in order for the RS model to provide a 
viable solution to the hierarchy problem. Quantum corrections to the Higgs potential in the 
RS model exhibit an even stronger divergence than in the SM [19]. However, as long as the 
Higgs sector is localized near the IR brane (i.e., in the vicinity oi t — 1), the effective cutoff 
■^uv(^) is of the order of several TeV, and hence the (big) hierarchy problem is solved. 

The question of how to introduce such a cutoff in practical one-loop calculations is far from 
trivial. However, at large loop momentum (of order several times Mkk) external momenta 
can be neglected, and hence there is a single 4D (euclidean) loop momentum pi = —p^ on 
which the cutoff should be imposed. We propose to associate a t-dependent cutoff with every 
vertex of a Feynman diagram. This can be thought of as modeling the effect of a form factor, 
which accounts for the impact of quantum gravity on energy scales above the effective Planck 
scale at that point. In general, the ti coordinates of the vertices are integrated over the 
entire bulk (e < < 1), and hence the cutoff values vary between the TeV scale and the 
fundamental Planck scale. This is indeed an important effect, which makes gauge-coupling 
unification possible in warped extra-dimension models [31, 32, 33, 34, 35]. The situation 
simplifies considerably for one-loop diagrams containing vertices with Higgs bosons, such as 
the one in Figure 1. Denoting the coordinate of the two gluons by ti and t2 and that of the 
Higgs boson by ts, the fact that ^3 ~ 1 ensures that the momentum cutoff on the 4D loop 
integral is 

Pe < min {Auv(ii), Auv(t2), Auv(*3)} = Axev ■ (67) 

The same mechanism guarantees that the hierarchy problem is solved in RS models. 

It is well known that evaluating the triangle diagram in Figure 1 in the 4D theory requires a 
consistent, gauge-invariant regulator such as dimensional rcgularization, the reason being that 
the integral is superficially UV divergent. If one were to evaluate it in D = 4 dimensions, then 
a momentum-independent term proportional to g'^'^ would appear, violating gauge invariance. 
In order to regularize the loop integral we thus evaluate it in D = 4 — 2e dimensions, with 
e > 0, so as to regularize UV divergences. We then find that the sum in (21) gets modified to 

E ^(^) , (68) 

where fi ~ AtcV is the rcgularization scale. For very large masses rriq^ ^ fi, the dimensional 
regulator gives rise to a suppression, which renders the sum over KK modes convergent even 
for arbitrary 0{1) Yukawa couphngs. The hmits N ^ 00 and 77 — > can therefore be taken 
without encountering any ambiguities. The contribution from the dangerous intermediate 
region from super-massive KK modes with m^^ ~ M^cak/v in Figure 2, which previously gave 
rise to an unsuppressed contribution of the form (65), now receives an extra suppression factor 
rf^, and vanishes when one takes the limit 77 ^ (at fixed e). The infinite sum then coincides 
with the result (31) up to harmless 0{e) corrections. 

In the dimensional rcgularization scheme, gauge invariance is manifest in the 4D theory. 
Since the dimensional regulator also regularizes the infinite KK sum, we arc guaranteed that 
the 5D theory remains gauge invariant, too. On the other hand, dimensional rcgularization 
is perhaps not the most intuitive way in which to introduce a UV cutoff. As an alternative, 
we will therefore rephrase the discussion in a rcgularization scheme based on using the hard 
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momentum cutoff given in (67). In order to ensure 4D gauge invariance in tliis case, we first 
extract two powers of tlie external gluon momenta by talcing appropriate derivatives, after 
which the remaining loop integral is superficially convergent. Introducing the UV cutoff on 
this integral, and neglecting the Higgs-boson mass compared with m^^, we obtain 



+ y 

For small masses nig^ <C Axev this reduces to the simple expression vg'j[^/mg^ used in the sum 
(21). For very large masses ruq^ ^ Axev, on the other hand, the UV cutoff gives rise to a 
strong suppression proportional to A^gy/m^^, implying that such heavy KK modes decouple 
rapidly. It follows that, due to physical reasons, the sum over KK modes in (21) is effectively 
truncated once the KK masses exceed the scale Atsv- In the RS model with a branc-localized 
Higgs sector, the scale M„cak/v at which the high- mass KK modes start to contribute a 
positive contribution to 'Eg^N.r]) with logarithmic growth is paramctrically much larger than 
the effective cutoff scale Axev- It is then appropriate to truncate the sum at a value ~ 
Atgy/Mkr corresponding to KK masses much smaller than M^eak/v- It follows that 

( AT 111 ax 2 \ 

where the truncation error has the form of a threshold correction, which is always present in 
effective-theory calculations. 

We can summarize the above discussion by emphasizing the subtle fact that, in order 
to obtain the correct answer for the gluon-gluon fusion cross section in the RS model, it is 
essential to employ a consistent UV regularization scheme when evaluating the loop integral, 
despite of the fact that this integral is convergent. When this is done, the convergence of 
the infinite sum appearing in (22) is improved in such a way that the order in which the two 
limits N ^ oo and 77 — > are taken becomes irrelevant. Specifically, in the two regularization 
schemes we have considered in our analysis, we find 



/ °° (7/2 



lim y 



3+6Af „ 



E 

n=l "^9" 



ri=0 



N^oo, ri^O ^ mg„ 
n=l ^" 



00 



E 

^ n=l 



(AT-max 2 



(71) 



In both schemes the infinite sums arc superficially convergent. The first line refers to the 
dimensional regularization scheme, while the second line corresponds to employing a hard 
momentum cutoff. It would be instructive to rederive these results by means of a properly 
UV-regularized 5D loop calculation, rather than by regularizing the sum over KK modes, as 
we have done above. In such a computation, in which there would be no reference to individual 
KK states, one should find that the correction to the hgg amplitude is of the form (71) once 
the regulator on the Higgs-boson profile is removed (r^ — > 0). For the simplified case of one 
generation, we have convinced ourselves that this is indeed the case. Details of this rather 
cumbersome analysis, as well as its extension to three generations, will be presented elsewhere. 
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The terms suppressed by a power of the UV cutoff, which appear on the right-hand side 
of (70), parameterize the difference between the asymptotic value j]^*-'^^'^^) of the infinite 
sum and the sum over a finite number of KK modes. More generally, such threshold terms 
may also arise from the yet unknown effects of the UV completion of the RS model. Prom 
a low-energy perspective, the only requirement on such a completion that is relevant to us is 
that it must cure the hierarchy problem, by taming loop momenta exceeding the fundamental 
scale of quantum gravity. As long as this is the case, the gluon-gluon fusion amplitude will 
also be regularized in the way discussed above. In the context of a specific UV completion, 
the threshold effects could be modeled at low energies by means of a brane-localized effective 
hG^^j^G^^'"' operator, suppressed by v/K\^^^. The additional factor reflects the high 

multiplicity of degrees of freedom in the low-energy effective theory. For quite generic reasons, 
the coefficient of this operator must contain the loop factor Q;s(/i)/(47r) factored out in (16), 
even in cases where the UV completion of the RS model is strongly coupled. The reason 
is that the on-shell external gluons couple proportional to their QCD charges, and that any 
new heavy state that couples to the Higgs boson must be color neutral, so it cannot have 
a tree-level coupling to gluons. Hence, a generic UV completion will indeed give rise to a 
threshold correction of the form shown in (70). The difference between the two sums (31) and 
(32), or the corresponding expressions (51) and (60), can however not be attributed to such a 
brane-localized threshold term. Rather, it is related to the question whether a physical cutoff 
and a sensible UV completion are present at all. 



8 Phenomenology 

We now study the implications of our results for Higgs-boson production in gluon-gluon fusion 
at the LHC. The master formula for the cross section has been given in (27), where the 
expressions for and Kg can be found in (17) and (28), respectively, while to an excellent 
approximation we can set Kg^ = 0. For the calculation of the Wilson coefficient we 
use our central result (60) and assume that the parametrically-suppressed threshold effects 
appearing in (70) can be neglected. Note that the trace over a function f{Xq) of the matrix 
defined in (11) is determined by the non-negative, real square roots yq'' of the eigenvalues of 
the hermitian matrices YgY^, i.e. 

Since in the RS model with anarchic 5D Yukawa couplings these matrices are structureless, 
the above result is proportional to the rank of the Yukawa matrices, which in the case of the 
minimal RS scenario is equal to the number of fermion generations. In other words, the KK 
towers of all six quarks give comparable contributions to the effective hgg vertex, irrespective 
of the mass of the corresponding SM fermion. We emphasize that this feature is not present in 
many other extra-dimensional extensions of the SM. For instance, in models based on universal 
extra dimensions the 5D Yukawa couplings are hierarchical, like in the SM, and hence the 
Higgs-boson couplings to light fermions and their KK excitations are strongly suppressed [36] . 
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More interestingly, in several gauge-Higgs unification models, in which the Higgs appears as a 
pseudo-Goldstone boson, one finds that the contribution to the hgg amplitude from the KK 
excitations of the SM quarks exactly cancels the dominant effect due to the corrections to the 
Yukawa couplings of the SM quarks, leaving only chirally-suppressed corrections, which are 
very small for all quarks other than the top quark [30].'^ Hence, in these new-physics scenarios 
only the top quark and its heavy partners contribute to the effective hgg couplings, while 
Higgs-boson production is independent of the masses and couplings of the KK excitations of 
hght SM quarks. 

By randomly generating a large set of 5D Yukawa matrices, which are required to satisfy 
^ I/max and to correctly reproduce the Wolfenstein parameters p and t] of the unitarity 
triangle, we have found that to a good approximation the average result for the trace appearing 
in (72) can be parametrized as^ 

Tr/(X,) ^ 3 (1 - 1.3 xLx) /(W) ; with f{x) = , (73) 

and Xmax = V Utaax/iV^MKK)- Likewise, the parameter Kt entering (28) can be approximated 
as 

Kt^l-l.lxl,,-eu, (74) 

where to an excellent approximation the contribution proportional to cancels in the sum 
C™(m/j) + Re (nt) ^(^t) appearing in the numerator of the expression for Kg in (28). Further- 
more, the very small effects due to Hh and ea can be neglected for all practical purposes. 

In the left plot in Figure 3, we show our results for the Higgs-boson production cross section 
in gluon-gluon fusion relative to the SM cross section, 

(^{gg h)sM 

as a function of the mass of the lightest KK gluon state, Mg(i) ^ 2.45 Mkk- We use the 
lightest KK gluon mass as a reference, because its value is a model-independent prediction of 
the RS scenario. The masses of the lightest KK fermions have very similar values but depend 
to some extent on the bulk mass parameters. Here and below we employ a Higgs-boson mass 
of 125 GeV. The solid red, green, and light blue lines show the approximate results obtained 
from (73) and (74) employing the values ?/max = 3, 1.5, and 0.5, respectively. The underlaid 
density plots indicate the distribution of the predictions for a large number of anarchic Yukawa 
matrices Yg, subject only to the constraint that their elements are smaller in magnitude than 
a given value ^max- The fact that the regions of highest density in the scatter plots nicely 
reproduce the results obtained directly from (73) and (74) shows that the requirement to 
reproduce the correct values of p and r] does not play an important role numerically. The 



^Thc authors of [37] argue, however, that this cancellation is model dependent and depends on the embed- 
ding of the fermions in the composite multiplets. 

*In [25], we preferred to use the associated Yukawa matrices Yg defined in (11) instead of the original 5D 
Yukawa matrices Yg. In practice, it turns out to be immaterial whether one considers or as random 
complex matrices. 
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Figure 3: Left: Predictions for the ratio Rh in the minimal RS model with bulk matter 
fields and an IR-Iocalized Higgs sector. The red, green, and blue density bands cor- 
respond to ?/max = 3, 1.5, and 0.5, respectively. The overlaid solid lines are obtained 
using the approximate parameterizations given in (73) and (74) for the same values of 
Z/max- Right: Contour plot for the ratio Rh obtained using the latter parameterizations. 



contour plot displayed on the right in the figure gives a two-dimensional representation for the 
cross section as a function of Mg(i) and |/max, obtained by employing again (73) and (74). 

We observe from Figure 3 that Rh is strictly below 1 and decreases (increases) with in- 
creasing i/max (KK scale). In other words, the minimal RS model with a brane-localized Higgs 
sector predicts a depletion of <j{gg — )■ h) relative to the SM. In the region where f i/max/^KK 
is a suitable expansion parameter, we obtain the approximate result 

Rh^ I -,^{U.2yL^ + 3.5), (76) 

where the constant term in parenthesis is due to the effect of k^. For not too small Yukawa 
couplings this observable is dominated by the effects of KK quark loops. Given the strong 
dependence of the ratio Rh on ?/max5 we find that for M^^^ ^ 3 TeV and Yukawa couplings close 
to the perturbativity bound |/max ~ 3 [18], the new-physics contributions to the Higgs-boson 
production cross section in gluon-gluon fusion can become so large that they completely can- 
cel the SM contribution. In fact, the sensitivity of Rh to the overall size of the 5D Yukawa 
couplings is even more pronounced than the one arising in the case of dipole-operator tran- 
sitions such as B ^ [38]. While the latter contributions also scale with y^^x; unlike Rh 
they are (at the one-loop level) insensitive to the multiplicity of states in the fermionic sector 
of the RS model under consideration. This feature underscores our assertion (made in the 
introduction) that precision measurements of the Higgs-boson properties furnish a superb tool 
for illuminating the quantum structure of electroweak interactions in RS scenarios. 



23 



In the minimal RS model considered here, constraints from electroweak precision observ- 
ables [39] and flavor physics [18] require that the lightest KK excitations of SM particles must 
have masses in the lOTeV range, which puts them outside of the reach for production at 
the LHC. Figure 3 shows that even in this case there can be significant virtual effects of KK 
particles on the Higgs-boson production cross section, provided that the 5D Yukawa couplings 
are not too small. The bounds from electroweak precision measurements, in particular, can 
be relaxed in several ways. For instance, "little RS models" , in which the size L of the extra 
dimension is reduced [40], would only have a minor impact on our analysis. Other extensions, 
such as models with a custodial SU{2)ii gauge symmetry in the bulk [9, 10], might however 
give rise to a rather different Higgs-boson phenomenology [14, 25] . 

9 Conclusions 

The announcements of the first direct hints for a Higgs-boson signal by the LHC and Tevatron 
experiments open up a new chapter in particle physics. Although the significance of the 
various measurements is not yet sufficient to preclude the possibility of statistical fluctuations 
accounting for the observed effects, hopes are high that, with ATLAS and CMS accumulating 
more data, a Higgs boson will be discovered (or excluded) by the end of this year. This 
discovery would not only mark the birth of the hierarchy problem, but it will also reshape 
some of the fundamental questions of our fleld. In particular, the focus of large parts of the 
LHC physics program will shift towards determining the Higgs-boson properties as accurately 
as possible, with the ultimate goal of probing possible deviations from the SM expectations. 
In close analogy to flavor physics, precision Higgs-boson physics represents a powerful way to 
investigate the dynamics of electroweak symmetry breaking at the quantum level. In fact, the 
couplings of the Higgs boson to a pair of gluons and photons vanish at tree level in the SM, 
but are induced by the exchange of virtual top quarks and W bosons at one-loop order. The 
effective hgg and hjj couplings thus offer a distinctive window to physics beyond the SM, 
where new heavy particles can propagate in the loops, thereby potentially affecting both the 
production cross section and the decay rates of the Higgs boson. 

The main goal of this article was to perform an analytic calculation of the impact of KK 
fermions on the production of a SM-like Higgs boson in the minimal RS scenario featuring an 
SU{2)l X U{1)y bulk gauge symmetry and a brane-localized Higgs sector. In particular, we 
have revisited the gluon-gluon fusion process, for which two independent calculations previ- 
ously found contradictory results. A significant suppression of the cross section was reported 
in [14], while in [15] an effect of similar magnitude but opposite sign was obtained. As we 
have shown, the discrepancy is not due to a simple computational mistake. On a technical 
level, it can be traced to the fact that, in order to make the overlap integrals of the wave func- 
tions of the brane-localized Higgs-boson with the bulk fermions mathematically well defined, 
a regularization of the Higgs profile is unavoidable in an intermediate step of the calculation. 
It is achieved by smearing out the profile over a finite width rj, thereby moving the Higgs 
boson slightly into the bulk. Given that the Yukawa couplings depend on rj, and that the 
calculations in [14, 15] both involve the summation over an infinite number N of KK levels, 
the question arises whether the result for the gg ^ h amplitude might depend on the order 
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in which the hmits 77 — )■ and iV — )■ 00 are taken. We have demonstrated that the order in 
which these hmits are performed indeed explains the aforementioned discrepancy. While the 
results of [14] are reproduced if one first takes the limit of the Higgs regulator to zero and 
then sums over the infinite tower of KK states, the findings of [15] correspond to the reversed 
order of taking the two limits. 

We have presented for the first time closed analytical expressions for both results in terms 
of the fundamental parameters of the RS model, valid to all orders in the ratio w/Mkk of the 
Higgs vacuum expectation value and the KK mass scale. We have then pointed out that the 
non-commutativity of the limits 77 — > and A?" — > 00 is due to a hidden UV sensitivity of the 
infinite KK sum, and illustrated this phenomenon by means of a simple toy model, which can 
be solved exactly. This computation highlights that for fixed 77 there exists a region, starting 
at around N ?a 0.1/?7, beyond which the KK-fermion wave functions resolve the width of the 
Higgs profile. As a result, the subtle cancellation of the Yukawa couplings of KK modes within 
each KK level, which is instrumental for the convergence of the KK sum, is spoiled. Thus 
logarithmic growth of the gg ^ h amplitude kicks in, which continues until the KK fermion 
profiles exhibit a large number of oscillations within the region where the Higgs profile is 
localized. Excitations with KK numbers higher than N pa 10 /r] quickly decouple, rendering 
the infinite KK sum ultimately convergent. 

The question which of the two calculations of the gluon-gluon fusion amplitude gives the 
correct physics result is resolved by introducing a proper UV regularization scheme. Specif- 
ically, we have studied the problem using dimensional regularization and using a hard mo- 
mentum cutoff, which may be identified with the inherent UV cutoff of RS models set by 
the warped Planck scale. In both cases, the regulator improves the convergence of the KK 
sum and removes the unphysical, previously unsuppressed contributions from super-heavy KK 
modes. The regularized sum over KK states yields the result of [14], up to possible threshold 
corrections suppressed by inverse powers of the effective UV cutoff. This shows that Higgs 
production is an UV-insensitive process in the RS model and can be calculated unambiguously 
once a proper regularization is employed. It would be worthwhile to reproduce our main result 
(60) for the induced hgg vertex by performing an actual 5D loop calculation rather than an 
infinite KK sum. This is left for future work. 

We have finally studied the numerical impact of our results for Higgs-boson production in 
gluon-gluon fusion, thereby extending our previous, more phenomenologically oriented analy- 
ses of Higgs physics in RS models [14, 25]. We have reemphasized the important point that, 
regardless of the mass of the corresponding SM fermion, in warped extra-dimensional models 
the loop-induced couplings of the Higgs field to gauge bosons receive similar contributions from 
the KK towers of each fermion state, and that these contributions scale with the square of the 
5D Yukawa couplings. Assuming a similar KK mass scale, RS scenarios hence predict much 
stronger effects than many other extra-dimensional extensions of the SM, such as universal 
extra dimensions and models in which the Higgs emerges as a pseudo-Goldstone boson. While 
we believe this fact to be true in general, it is not unlikely that the structure of the corrections 
depends on the precise realization of the scalar sector. Exploring the associated model depen- 
dence in more detail is left for future work. In order to probe the large effects predicted by 
our calculations, it would be particularly useful to measure the cross sections for Higgs-boson 
production in both the gluon-gluon and vector-boson fusion channels. A combination of theses 
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two measurements would allow for a clean extraction of the gg ^ h amplitude, and hence 
represents a unique way to study the Yukawa sector of RS theories. 
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A Special case of vanishing bulk masses 

In the case where all bulk mass parameters ca arc set to zero, and where in addition we set 
e = (which in this case is unproblematic), the EOMs (7) can be solved exactly if we work with 
a sufficiently simple regularized Higgs profile. This provides a nice test case for our general 
results and conclusions. In order to obtain analytic expressions for the wave functions of the 
individual KK modes, we consider the Higgs profile (61) and furthermore assume that there is 
a single fermion generation with real and positive Yukawa coupling Yq. The generalized mass 
matrix in (8) is then given by 



vY„ 1/0 1 



Mg{t) = <^ v^Mkk V\l 



= z 



1 

1 



0; 



for 1 — T] < t < 1 , 
for e < t < 1 — rj . 



(Al) 



The wave functions Q^a\^) ^'^'^ 2-component objects. 

In the solution of the EOMs we must distinguish two cases. For Xn < z, we find 



Qin) 



it) 



( cos(x„t) \ 
T sin(x„t) 
sin(x„i) 
V ± cos(x„i) / 



t < l-?7. 



(A2) 



and 



Ql\t) 



/ ri cosh [■\/z'^ - 






sinh [s/ z'^ 




-*)] 


r2 sinh \\/ - 






\ ±ri cosh \\J z^ 


-4(1 


-t)]) 



t > 1-77, 



(A3) 
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where the coefficients are given by 



ri 



cos [,/•„( 1 - //)] 



sill 



cosh [y^ 
For Xn> z,we find instead 



r2 



['■nil 



sinh [^z'^-xIt]] 



(A4) 



Qin) 



Nr. 



I ri cos [a/x^ - z^(1- t)] ^ 
=Fr2 sin [^/xl - z"^ (1 - t)] 
r2 sin [y/a;2 - (1 - t)] 

V in cos [y^2;3^(i_^)]y 



t>l-r;, 



(A5) 



where 



cos [xn{l - v)] 



sin [xn{l - v)] 



= 2 2^ ' = ^^^^===^ . (A6) 

cos [ V - -z^ 77J sm [ yx^ - 2;^ 77J 

The eigenvalues x„ are determined by matching the profile functions at t = 1 — r^. This 
leads to the conditions shown in (62). Evaluating the normalization constraint, wc obtain the 
expression for given in the second line of (63). Using these results, it is not difficult to 
calculate the relevant overlap integrals with the Higgs profile. In this way we arrive at the 
result for the sum T,q{N, rj) shown in the first line of (63). 
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